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Abstract. Conformally flat tilted Bianchi type V cosmological models in presence 
of a bulk viscous fluid and heat flow are investigated. The coefficient of bulk viscosity 
is assumed to be a power function of mass density. Some physical and geometric 
aspects of the models are also discussed. 
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1. Introduction 

A considerable interest has been shown to the study of physical prop- 
erties of spacetimes which are conformal to certain well known gravi- 
tational fields. The general theory of relativity is believed by a number 
of unknown functions - the ten components of g 13 . Hence there is a 
little hope of finding physically interesting results without making re- 
duction in their number. In conformally flat spacetime the number of 
unknown functions is reduced to one. The conformally flat metrices 
are of particular interest in view of their degeneracy in the context of 
Petrov classification. A number of conformally flat physically significant 
spacetimes are known like Schwarzschild interior solution and Lemaitre 
cosmological universe. 



The study of Bianchi type V cosmological models create more inter- 
est as these models contain isotropic special cases and permit arbitrarily 
small anisotropy levels at any instant of cosmic time. This property 
makes them suitable as model of our universe. Also Bianchi type V 
models are more complicated than the simplest Bianchi type models 
e.g. the Einstein tensor has off diagonal terms so that it is more natural 
to include tilt and heat conduction. Spacetimes of Bianchi type I, V and 
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IX universes are the generalizations of FRW models and it will be inter- 
esting to construct cosmological models of these types which are of class 
one. Roy and Prasad (1994) have investigated Bianchi type V universes 
which are locally rotationally symmetric and are of embedding class one 
filled with perfect fluid with heat conduction and radiation. Bianchi 
type V cosmological models have been studied by other researchers 
(Farnsworth, 1967; Maarteins and Nel, 1978; Wainwright it al. , 1979; 
Collins, 1974; Meena and Bali, 2002). 

The general dynamics of tilted models have been studied by King 
and Ellis (1973) and Ellis and King (1974). The cosmological models 
with heat flow have been also studied by Coley and Tupper (1983. 
1984); Roy and Banerjee (1988). Ellis and Baldwin (1984) have shown 
that we are likely to be living in a tilted universe and they have in- 
dicated how we may detect it. Beesham (1986) derived tilted Bianchi 
type V cosmological models in the scale-covariant theory. A tilted cold 
dark matter cosmological scenario has been discussed by Cen, Nickolay, 
Kofman and Ostriker (1992). Recently Bali and Meena (2002) have 
investigated two tilted cosmological models filled with disordered radi- 
ation of perfect fluid and heat flow. Tilted Bianchi type I cosmological 
model for perfect fluid distribution in presence of magnetic field is inves- 
tigated by Bali and Sharma (2003). Tilted Bianchi type I cosmological 
models filled with disordered radiation in presence of a bulk viscous 
fluid and heat flow are obtained by Pradhan and Rai (2003). 

Most cosmological models assume that the matter in the universe 
can be described by 'dust'(a pressureless distribution) or at best a 
perfect fluid. Nevertheless, there is good reason to believe that - at least 
at the early stages of the universe - viscous effects do play a role (Israel 
and Vardalas, 1970; Klimek, 1971; Weinberge, 1971). For example, the 
existence of the bulk viscosity is equivalent to slow process of restoring 
equilibrium states (Landau and Lifshitz, 1962). The observed physical 
phenomena such as the large entropy per baryon and remarkable de- 
gree of isotropy of the cosmic microwave background radiation suggest 
analysis of dissipative effects in cosmology. Bulk viscosity is associated 
with the GUT phase transition and string creation. Thus, we should 
consider the presence of a material distribution other than a perfect 
fluid to have realistic cosmological models (see Gr0n ) for a review 
on cosmological models with bulk viscosity). The model studied by 
Murphy (1973) possessed an interesting feature in that the big bang 
type of singularity of infinite spacetime curvature does not occur to be a 
finite past. However, the relationship assumed by Murphy between the 
viscosity coefficient and the matter density is not acceptable at large 
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density. The effect of bulk viscosity on the cosmological evolution has 
been investigated by a number of authors in the framework of general 
theory of relativity .(Pavon, 1991; Padmanabhan and Chitre, 1987; 
Johri and Sudarshan, 1988; Maartens, 1995; Zimdahl, 1996; Santos et 
ai, 1985; Pradhan, Sarayakar and Beesham, 1997; Kalyani and Singh 
1997; Singh, Beesham and Mbokazi, 1998; Pradhan et al, 2001, 2002, 
2003). This motivates to study cosmological bulk viscous fluid model. 
Recently Meena and Bali (2002) have investigated two conformally flat 
tilted Bianchi type V cosmological models filled with a perfect fluid and 
heat conduction. In this paper, we propose to find tilted Bianchi type 
V cosmological models in presence of a bulk viscous fluid and heat flow. 



2. The metric and field equations 

We consider the Bianchi type V metric in the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 e 2x (dy 2 + dz 2 ) , (1) 
where A, B are function of t only 

The Einstein's field equations (in gravitational units c = 1, G = 1) 
read as 

i^-^ = -87rT/, (2) 

where R\ is the Ricci tensor; R = g^Rij is the Ricci scalar; and T- is 
the stress energy-tensor in the presence of bulk stress given by 

Ti = (p + p)viV j + pgi + q iV j + Vi q> , (3) 

and 

P = P- £4- ( 4 ) 

Here p, p, p and £ are the energy density, isotropic pressure, effective 
pressure and bulk viscous coefficient respectively and Vi is the flow 
vector satisfying the relations 

9ijv i v j = -1, (5) 

m j > o, (6) 

qy = 0, (7) 

where % is the heat conduction vector orthogonal to v\. The fluid flow 
vector has the components ( S1 °^ A , 0, 0, cosh A) satisfying Eq. (5) and A 
is the tilt angle. 
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The Einstein's field equations (2) for the line element (1) has been set 
up as 

-8vr[(p+p)sinh 2 A + p + 2^ 1 sinhA] = (j^j - (8) 

-8vr[-(p+p)cosh 2 A+p-2^ 1 sinhA] = -^+^J (10) 

— 87r[(p+p)y4sinh A cosh A+^4 2 gi(cosh A+sinhA tanhA)] = — — , 

(11) 

where the suffix 4 at the symbols A, B denotes ordinary differentiation 
with respect to t. 



3. Solution of the field equations 



Equations (8) - (12) are four independent equations in seven unknowns 
A, B, p, p, £, q and A. For the complete determinacy of the system, we 
need three extra conditions. 

First we assume that the spacetime is conformally flat which leads to 



C2323 — - 
and secondly, we assume 



^44 -B44 Aa^B^ B^ 

~A B AB~ + 52 







A = B n , 

where n is any real number. Eqs. (12) and (13) lead to 



-B44 . ,-B 2 
-B- + {n - 1) m 

From Equations (8), (10) and (13), we have 



0. 



-4vr [(p + p) cosh 2A + 4B n qi sinh A] = ^ 

B 



B\ 1 

71 -=-77 + 



B 2 B 2r, 



and 



. . -B44 , . B A 2 



Equations (9), (16) and (12) lead to 

/ 2 , 1^4 ^44 

{n-n + l)—- n — 



B 2n 



4ir(p+p). 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 
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Equations (11) and (13) lead to 

™ , > 2(n - l)_B 4 tanh2A . . . , 
-167r^iS n smhA = + 4n(p + p) sinh 2A tanh2A. 

(18) 

From Eqs. (15) and (18), we obtain 

-B 44 nB\ 1 4ir(p + p) 2(n-l)tanh2A 

~B B~ + ~B 2 ~™~ cosh2A + W l ' ( ' 

Equations (17) and (19) lead to 

Bu nB\ 1 _ nB 44 (n 2 - n - \)B\ g2 „ 2(n - 1) £ 4 tanh2A 

5 5 S 2 ™ B B 2 + sec + Bn+1 

(20) 

Equation (14) can be rewritten as 

which on integration leads to 

B = (at + (3)n. (22) 
where a, (3 are constants of integration. Hence we obtain 

A 2 = n 2 {at + P) 2 , (23) 

B 2 = nn (at + /3)n. (24) 
Hence the geometry of the spacetime (1) reduces to the form 

ds 2 = -dt 2 + n 2 (at + I3fdx 2 + [n(at + (5)]^e 2x {dy 2 + dz 2 ). (25) 

After the suitable transformation of coordinates, the metric (25) takes 
the form 

ds 2 = + n 2 T 2 dX 2 + nlTle 2X {dY 2 + dZ 2 ). (26) 

The effective pressure and density of the model (26) are given by 

8vrp = 8^-^) = -^^, (27) 

torn 3 (» 2 ~ 1) 

87T P = n 2 T 2 ' ( 28 ) 
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where is the scalar of expansion calculated for the flow vector v l and 
given is as 

e= 2K + (n + 2)ka 
nT 

The tilt angle A is given by 

cosh 2 A = fc 2 , (30) 

sinh 2 A = K 2 , (31) 
where k and K are constants given by 

k 2 = (na 2 -!) 2 

(a 2 - l)(2n 2 a 2 - na 2 - 1) + 2a 2 (n - 1) ^(n 2 - n)(l - a 2 ) ' 

(32) 

r2 _ (n - l)(2raa 2 - na 4 - a 2 ) - 2a 2 (n - 1) ^{n 2 - n){\ - a 2 ) 

" (a 2 - l)(2n 2 a 2 - na 2 - 1) + 2a 2 (n - l)^(n 2 - n)(l - o?)' 

(33) 

If we put £ = in (27), we get the solutions as obtained by Meena and 
Bali (2002). 

Thus, given £(t) we can solve the system for the physical quantities. 
Therefore to apply the third condition, let us assume the following 
adhoc law (Maartens, 1995; Zimdahl, 1996) 

m = iop m (34) 

where £o and m are real constants. If m = 1, Eq. (34) may correspond 
to a radiative fluid (Weinberg, 1972), whereas m = \ may correspond to 
a string-dominated universe. However, more realistic models (Santos, 
1985) are based on lying the regime < m < \. 

3.1. Model I: (£ = £ ) 

When m = 0, Equation (34) reduces to £ = £o = constant and hence 
Equation (27) with the use of (29) leads to 

_ g {2K + (n + 2)fc«} (a 2 - 1) 
P ~ nT SrfT 2 ' (35j 



3.2. Model II: (£ = £ p) 

When m = 1, Equation (34) reduces to £ = £oP and hence Equation 
(27) with the use of (29) leads to 



1 

P 



8irn 2 T 2 



1 _ a 2 3tjo(a 2 -l){2K + (n + 2)ka} 
a nT 



(36) 
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It is observed from Equations (28), (35) and (36) that p and p vary as 
7^. The models are singular at T = and as they evolve, the pressure 
and density decrease. 



4. Some Physical and Geometric Properties of the Models 

The weak and strong energy conditions, we have, in Model I 



_ (a 2 - 1) £q{2K + (n + 2)ka} 
^ ^ 4:irn 2 T 2 nT 

(a 2 - 1) ^{2K +{n + 2)ka} 



P-P 



p + Sp 



2irn 2 T 2 nT 

3£ {2K + (n + 2)ka} 



nT 



p-3p 



(3a 2 - 1) 3£ {2K + (n + 2)ka} 
Aim 2 T 2 nT 



In Model II, we have 
1 



P + P 



P-P 



87rn 2 T 2 
1 

8vm 2 T 2 



2{a A - 1) + 



4(a 2 - 1) - 



3^ (a 2 - 1){2K + (n + 2)/ca} 

3^ (a 2 - l){2if + (n + 2)ka} 
nT 



p + 3p 



9£ (" 2 - 1){2K + (n + 2)fea} 



p-3p = 



8vrn 2 T 2 



2(a 2 - 1) + 



8vm 3 T 3 
3g (a 2 - 1){2K + (n + 2)fca} 



(37) 

(38) 
(39) 
(40) 

(41) 

(42) 

(43) 
(44) 



The reality conditions p > 0, p > and p — 3p > impose further 
restrictions on both of these models. 

The flow vector v 1 and heat conduction vector q l for the models (26) 
are obtained by Meena and Bali (2002) 



K 



v 4 = k, 



k{(a 2 - l)kK+(n-l)a) 
4irn 3 T 3 (k 2 + K 2 ) ' 



(45) 
(46) 
(47) 
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4 _ K{{a 2 -l)kK + {n-l)a} 
q 47rn 3 T3(fc 2 + i^ 2 ) ' { ' 

The rate of expansion Hi in the direction of X, Y, Z-axe are given by 

Hi = ^, (49) 

H2 = H 3 = -£L (50) 

The non- vanishing components of shear tensor (cr^ ) and rotation tensor 
(wjj) are obtained as 

a u = -nk 2 T[(n- l)ka - K], (51) 
3 

V22 = a 33 = [(1 - n)ka - 2K], (52) 

2K 2 

a A4 = —[(n-l)ka-K], (53) 
Snl 

cr±4 = —[2(1 — n)k 2 + 2kK — 3n], (54) 

wi4 = naK. (55) 

The models, in general, represent shearing and rotating universes. The 
models start expanding with a Big bang at T = and the expansion in 
the models decreases as time increases and the expansion in the models 
stops at T = oo and a = — ^+2)k • Both density and pressure in the 
models become zero at T = oo. For a = 1, n = 1, we observe that heat 
conduction vector q 1 = q 4 = 0. When T — > oo, v 1 = 0, w 4 = constant, 
g 1 = g 4 = 0. Since lim^oo f 7^ 0, the models do not approach isotropy 
for large values of T. There is a real physical singularity in the model 
at T = 0. 



5. Particular Models 

If we set n = 2, then the geometry of the spacetime (26) reduces to the 
form 

dT 2 

ds 2 = — + TdX 2 + Te x (dY 2 + dZ 2 ). (56) 

The effective pressure and density for the model (56) are given by 

87rp = 87r(p-^) = -^H, (57) 
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where 6 is the scalar expansion obtained as 

e = ^L±^. (59) 

The tilt angle A is given by 

cosh 2 X = kf, (60) 

sinh 2 A = Kl, (61) 
where k\ and K\ are constants given by 

1 6a 4 - 7a 2 + 1 + 2aV2(l - a 2 ) ' V ; 



a 2 

if i 2 _ 



(3 - 2a 2 - 2^2(1 -a 2 ) 



L1 6a 4 -7a 2 + l + 2a 2 V2(l-a 2 )' (63) 
Thus, for given £(t) one can solve the system for the physical quantities. 

5.1. Model I: (£ = £ ) 

When m = 0, Equation (34) reduces to £ = £o = constant and hence 
Equation (57) with the use of (59) leads to 

Z {K 1 + 2k 1 a} (« 2 -l) 
P T 32vrT 2 " { ' 

5.2. Model II: (f = £ p) 

When m = 1, Equation (34) reduces to £ = £o/0 and hence Equation 
(57) with the use of (59) leads to 



1 

P ~ 32vrT 2 



. 2 , 3^(a 2 - l){#i + 2k ia } 
I — a H — 



(65) 



It is observed from Equations (58), (64) and (65) that p and p vary as 
The models are singular at T = and as they evolve, the pressure 
and density decrease. 
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6. Some Physical and Geometric Properties of Particular 

Models 



The weak and strong energy conditions, we have, in Model I 

* 2 -l) Z {K 1 + 2k 1 a} 



P + P 16vrT 2 ' T 

_ {a 2 - 1) _ gojifi + 2fcia} 

P P ~ 8vrT 2 T 

3£o {K! + 2k ia } 
p + 3p = , 

_ (3a 2 - 1) 3g {Ki + 2k ia } 
9 6P ~ 16ttT 2 r 



In Model II, we have 

P + P= V^ 
l 



P-P 



32vrT 2 



2 , , 3Co(a 2 -l){K 1 +2fc 1 q} 
If' 2 1) , 



2(a 2 -l)+ r 

.2 ^ 3g (a 2 -l){if 1 +2fc 1 a} 



_ 9e (a 2 -l){Jri + 2fcia} 



p-3p = 



32vrn 2 T 2 



2(a 2 - 1) + 



3^ (a 2 -l){K 1 +2fci«} 



(66) 

(67) 
(68) 
(69) 

(70) 

(71) 

(72) 
(73) 



The reality conditions p > 0, p > and p — 3p > impose further 
restrictions on both of these models. 

The flow vector v l and heat conduction vector q l for the models (56) 



are obtained as 



Q = - 



4 _ _ 



i Ki 
V =2T' 

v 4 = k\, 

ki{(a 2 - l)kiKi + a} 
327rT 3 (fc 2 + K 2 ) ' 

i^i{(a 2 - 1)^1^1 + a} 



(74) 
(75) 
(76) 

(77) 



327rT 3 (fc 2 + K 2 ) 

The rate of expansion Hi in the direction of X, Y, Z-axe are given by 

(78) 



#1 = -, 

± rp ' 
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H 2 = H 3 = (79) 

The non- vanishing components of shear tensor (<7y) and rotation tensor 
(uJij) are obtained as 

* 11 = lk%T[k 1 a-K 1 ], (80) 

IX 

<722=<733 = -^g-[fcia + 2#l], (81) 

K 2 

<744 =g^[fcia-^i], (82) 

<7i4 = [ 2 *t - 2fei ifi + 6], (83) 

wi4 = 2aKi. (84) 

The models in general represent shearing and rotating universes. The 
models start expanding with a Big bang at T = and the expansion 
in the models decreases as time increases and the expansion in the 
models stops at T = oo and a = — Both density and pressure in 
the models become zero at T = oo. For expansion and rotation, we 
have a ^ 1, K± ^ 1 when T — > oo we observe that heat conduction 
vector q l = q A = 0, v 1 = 0, v A = K\. Since lim^oo f / 0, the models 
do not approach isotropy for large values of T. There is a real physical 
singularity in the model at T = 0. When k± = 1 then A = 0. Thus, the 
tilted cosmological models lead to non-tilted one for k± = 1. 



7. Conclusions 

We have obtained a new class of conformally flat tilted Bianchi type V 
magnetized cosmological models with a bulk viscous fluid as the source 
of matter. Generally, the models are expanding, shearing and rotating. 
In all these models, we observe that they do not approach isotropy for 
large values of time T in the presence of magnetic field. 
The coefficient of bulk viscosity is assumed to be a power function of 
mass density. The effect of bulk viscosity is to introduce a change in 
the perfect fluid model. We also observe here that the conclusion of 
Murphy (1973) about the absence of a Big bang type of singularity in 
the finite past in models with bulk viscous fluid is, in general, not true. 
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